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Abstract 

Using 9.0 fb -1 of integrated luminosity in e + e _ collisions near the T(4S) mass collected with the 
CLEO II. V detector we report the first observation of the decay D° — > K^nir . We measure the 

ratio of branching fractions, ^^jol^o^o ^ = 0-46 ± 0.07 ± 0.06. We perform a Dalitz analysis of 

155 selected D° — > Kgnir candidates and find leading contributions from ao(980)K<j and K*(892)n 
intermediate states. 



*On leave of absence from University of Chicago. 
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A large fraction of the known D meson decay rate is in three-body hadronic decays to 
the pseudoscalar particles K and tt. These decays dominantly proceed through quasi-two- 
body intermediate states with a rich set of resonances. The dynamics of three body decays 
can be studied using the Dalitz technique Interest in the decay D° — > K s r]n stems 
from comparing the results of the Dalitz plot analyses of the decay D° — > K s tt + 7t~ studied 
by ARGUS and CLEO0 with the decay D° -> K S K + K~ studied by AUGUS T and 
BaBar[5|. The contribution of the /o(980) observed in the former case is not enough to 
explain the ~60% fraction observed in the latter decay. Additional scalar contribution from 
a o (9S0)K s can be expected in D° -> K S K + K', but is difficult to separate from f (9S0)K s 
in the Dalitz plot. The a (980)K s intermediate state can also be observed in the favored 
a (980) — > i]tt° decay mode which would give rise to the D° — > K s r]ir final state. The decay 
D° — > Kg7]7T° or any other D° modes with ao (980) in the intermediate state have not yet 
been observed. There is little information on D° decay modes with rj in the final state; only 
an upper limit BR(D° -> r)X) < 13% @ C.L. = 90% has been measured. Note that 
Kgi]7T is a CP eigenstate. A large sample with a good signal to noise ratio in this mode 
can be used for studies of CP violation in D° and D° decays. 

The data sample used in this analysis was produced by the Cornell Electron Storage 
Ring (CESR) and collected with the general purpose CLEO II. V (7( detector. Our analysis 
is based on 9.0 fb _1 of integrated luminosity of e + e~ collisions at ^fs ~ 10 GeV above and 
below BB production threshold. Charmed particles can be produced both in the process 
e + e~ — ► cc and in B meson decays. To suppress events with low momentum D 0, s from B 
decays, which have higher multiplicity and higher combinatorial backgrounds, we use the 
decay D* + — ► D tt + (charge conjugation is implied throughout this letter) as a tag and 
require that the D* + momentum exceeds 2.8 GeV/c. The decay D° — ► K s r]7T° is observed 
in the most probable mode of the final state, K s — > 7r + 7r~, i] — ► 77, 7r° — > 77. 

Charged tracks are required to be well measured in the tracking detectors. Candidate 
Kg-s are reconstructed from pairs of oppositely charged tracks assumed to be pions. The 
candidate Kg trajectory is required to be consistent with production in the interaction 
region, while its vertex should be significantly (> lOer) isolated from this region. We select 
Kg candidates if the reconstructed mass, m n + n -, is within 10 MeV/c 2 of the nominal K s 
mass p. On average, Kg-s in this selection have a mass resolution of a K o = 3.7±0.2 MeV/c 2 . 

We form 7r° and rj candidates from pairs of neutral showers in the CLEO Csl calorimeter. 
They are required to be consistent with electromagnetic showers, have an energy deposition 
above 30 MeV and be in the central, barrel region of the detector. For 7r° candidates 
we require the invariant mass, m 77 , of the photon pair to be within 18 MeV/c 2 of the 
nominal 7r° mass 0. The average detector resolution of 7r° — > 77 invariant mass is <J n o = 
6.1 ± 1.2 MeV/c 2 . Similarly, 77 candidates are required to have a two photon invariant 
mass within 40 MeV/c 2 of the nominal 77 mass [6] at the average detector resolution of 
a v = 12.6 ± 1.0 MeV/c 2 . 

We kinematically fit Kg, tt° and 77 candidates and constrain their masses to nominal 
values. This procedure improves the D° mass resolution by a factor of two for D° — > 
K s i]tt° decays. We reconstruct D° —>■ K s i]tt candidates by combining the Kg, tt° and rj 
candidates in the event. To eliminate the significant combinatoric background, we select the 
combination with the smallest 

Xm = ( m " 1 ~'~ m7] j + ^ m T"^ ^0 j _|_ — £.j ^ where a \\ the invariant masses are taken 
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FIG. 1: Distribution of the energy release, Q, in the decay D* 
difference of D — > Kg-qir candidates (b). 



D°tt + (a), and the mass 



before the mass constraint of the kinematic fit. 

The D° candidate is combined with tt + tracks to form the tagging decay D* + — > D°tt + . A 
significant D° signal is observed both in the energy release, Q = m(Kgr]7i°'n- + ) — m(Kgi]7T°) — 
m w +, and in the D° mass difference Am = m(Kgi]7i°) — m^o shown in Fig. ^ The Q 
distribution, shown in Fig. represents raw Q vs Am events in ~ 3cr signal Am band 
indicated by arrows in Fig. and vice versa. 

We estimate a signal yield of 155±22 events from a fit with a single Gaussian for the sig- 
nal plus a linear background to the mass spectrum of Fig. ^b). The GEANT-based Monte 
Carlo simulation (if of the CLEO II. V detector response is used to estimate the efficiency 
s(D° -> K° s r)7c°) = (1.15 ± 0.05 ± 0.12 ± 0.01)%, where the uncertainties are statistical, 
systematic and from the uncertainties on the Kg — > 7r + 7r~, rj 77, ir° 77 branching 
fractions, respectively. The systematic uncertainty includes the track reconstruction effi- 
ciency (2%/track), tc° and r] selection (5% each), and the background subtraction in D° 
mass spectrum (7.2%). The first two uncertainties absorb a variation of efficiency between 
the phase space and the resonant event production mechanism. The background subtraction 
error is estimated from variation in the signal yield when we change the fit function includ- 
ing a single versus double Gaussian for the signal and background described with a linear 
function, taken from the D° mass spectrum sidebands, or taken from the Q distribution 
sidebands. 

To measure the branching fraction we normalize to the total number of D 0, s pro- 



K° s 7i° 



decay with known rate, 
We use the same selection 



duced in the decay D* + —> D°n + . We use the D° 

BR(D° -> K° s n°) = lBR(D° -> K°tt ) = (1.14 ± 0.11)% 

as D° — > K^rjir , but without the rj, and find a very clean D° — > K^ir signal with yield of 

1105±54 events and an efficiency e(D° -> K° s 7r°) = (3.76±0.18±0.26±0.02)%. We find the 

►ifV°) _ 
^=0 — " 



ratio of branching fractions to be BR( - D 



BR(D°^K°rrn-°) 



0.46±0.07±0.06±0.003 



BR(D°^Wn°) BR(D°^ S „ 

0.46 ± 0.09, where the errors are statistical, systematic, and 77 — > 77 branching fraction 
uncertainties, respectively. Using the known D° — > K^it branching fraction, we find 

BR(D° -> Tfr/TT ) = (1.05±0.16±0.14±0.10)%, where the last error is associated with the 
uncertainty on the D° — > Kgir branching fraction. Many systematic uncertainties cancel in 
the ratio measurement. 

The selected sample, although small, is clean enough to search for possible intermediate 
states using the Dalitz technique jjj]. We tighten the mass difference selection criteria to 
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FIG. 2: Dalitz plot of D° — ► K^ijir (a), and the map for the adaptive binning (b). 



two standard deviations (|Am| < 25 MeV/c 2 and |AQ| < 1.2 MeV/c 2 ) in order to increase 
signal to background ratio. We select for Dalitz analysis 155 events (accidentally the same 
number of events that we find for measurement of the branching ratio) shown in Fig. |2Ja) as 
m 2 (r]7r ) versus m 2 (KgiT ). The same selection criteria were applied to measure the efficiency 
across the Dalitz plot with a simulation of D° —>■ Kgrj7r° decaying uniformly in its allowed 
phase space. The shape of the small background is taken from the data sample of 171 
events in a Q sideband, 10 < Q < 25 MeV/c 2 , and an extended range of invariant mass, 
|Am| < 100 MeV/c 2 . Both the efficiency and the background are nearly uniform across 
the Dalitz plot and we parameterize them separately with a two-dimensional polynomial of 
third degree obtained from the dedicated fit. 

The Dalitz plot, Fig. Ufa), shows a significant contribution from ao(980)Kg interfering 
with other resonances, as evidenced by the deficit of event density in the center of the plot, 
and by the shift to the left (right) of the ao(980) band on the top (bottom) of the plot. There 
is an indication of a K * (892)77 contribution as there is an enhancement in the expected region 
in the m 2 (Kgir ) projection, shown in Fig. |3J^a). The visible mass peak is shifted lower than 
would be expected given the if* (892) mass indicating interference of K*(892)r] with other 
intermediate states. 

To extract information from the Dalitz plot we apply the technique developed in our 
previous analyses [3] which uses an unbinned maximum likelihood fit and an "isobar 
model" to measure matrix element amplitudes. An isobar model approximates the matrix 
elements as 

M = a NR e^N R +j: R a R e^nAj({K s ,7 1 ,n }\R), 

a coherent sum of non-resonant (NR) and resonance (R) terms, each multiplied by its own 
complex factor. The complex factor is parametrized by a real amplitude a R and a phase 
ip R , which are extracted from the fit. The amplitude, Aj(ABC\R), is defined for the decay 
chain D° — > RC — > ABC with an intermediate resonance R represented by the Breit-Wigner 
function with spin J dependent factor. The overall amplitude normalization and complex 



phase are arbitrary, and are chosen such that a, 



do 



0. 



The mass m dependent width of the ao(980) is parameterized using the method of 
Flatte while the partial width is proportional to the phase space factor p = 2p/m 
instead of the decay momentum p, 



™r ao(980) (m) 

for the coupling constants 
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We assume an isospin symmetry 
9 2 aoK °K+/ 2 - In 
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m(a (980)) = 999 ± 5 MeV/c 2 
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O = 11-1 ± 1-0 GeV 2 , g 2 aoK+K -/g 2 aovn0 = 0.58 ± 0.09. 

The event density of the Dalitz plot is fit to the efficiency corrected matrix element 
squared and the background polynomial which is added incoherently [^] to the signal. The 
relative signal fraction, 0.867±0.027, is estimated from the Am spectrum of the data sample. 
In all Dalitz fits the signal fraction is a parameter of the fit constrained to this estimate. 

With our sample we find the most reliable goodness of fit estimator to be a x 2 -like 
parameter for Poisson statistics: x 2 = YliLi ( ~ ni ~ u ^ , |6[ where rij, ^ are the number of events 
and its mean expectation in the i-th bin, and N is a total number of bins. We split the Dalitz 
plot into 10 x 10 equal bins. In order to provide sufficient statistics for a mean expectation, 
it was necessary to join some bins using so-called "adaptive binning" and requiring vi (or 
rii)>5 in each bin. The 24 bins found with adaptive binning are shown in Fig. Efb). We 
have tested this goodness of fit parameter in simplified Monte-Carlo based simulations of 
our data and find that it gives a uniform probability for statistically distributed data with 
P(x 2 /Nd.o.f.) in the range [0,1]. The simulation of different models shows that we are only 
sensitive to contributions to the Dalitz plot that are greater than 20% of the total rate, 
and thus our goal is to find a consistent description of the observed event density using a 
minimal set of dominant modes. 

From previous observations we expect r/7r to have contributions from intermediate 
states including a (980), a2(1320), and ao(1450). Similarly, K^tt should have contributions 
from K*(892), if* (1410), #5(1430), K*(U30), and 7^(1680). A possible low-mass Kn- 
scalar state or dynamical structure, n, which is not included in [f| but is widely discussed in 
recent publications [12], could also contribute. There is no obvious contribution from KgT] 
in this mass range. We start with a minimal set of resonances and recognize an additional 
resonance as contributing if the fit probability improves, the amplitude is at least three 
standard deviations from zero, and the error on the phase is less than 30°. 

We find that a model including only a (980)Kg and K*(892)r] contributions gives a low 
probability of 0.8% and is an unlikely explanation of our data. Models with a single resonance 
are even worse with probabilities of less than 10~ 6 . Good consistency with our data can 
be achieved with models including two main intermediate states, ao(980)Kg, K*(892)i], and 
additional mode(s). We find four additional modes giving a fit probability >1%: (i) a non- 
resonant fraction; (ii) ^(1430)^; (iii) _K"q(1430)?7 and 02(1320) X° (fit projections are shown 
in Fig. I2J); and (iv) a k with parameters taken from [l2j. We do not find any significant 
contribution or fit quality improvement by adding other resonances. For these four models 
Table summarizes the amplitude and phase we extract from the fit for the .fT* (892)7/ mode, 
fixing the amplitude and phase for the ao(980)Kg mode to be one and zero respectively. 
Our sample is too small to allow us to choose one model among these four. In the last row 
of Table U we present averaged results and their variation due to our inability to choose a 
single decay model that describes our data adequately. 

When the amplitudes and phases are extracted from the fit we derive the fit fraction (FF) 
for each contribution. The fit fraction is defined for each resonance as its matrix element 
amplitude squared (rate) integrated over the allowed phase space divided by the total matrix 
element amplitude squared integrated over the same phase space. In general the sum of the 
fit fractions does not have to equal one due to interference among the contributions. A 
statistical uncertainty on the fit fraction is computed from the fit covariance matrix using 
Monte Carlo methods as described in @. Table |U gives the fit fractions for the ao(980)Kg 
and K*{892)rj modes, and the fit fraction for the additional mode(s), their averaged values, 
and estimated variations due to the choice of decay model. 
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TABLE I: Results for four models of the additional contribution beyond ao(980)Kg and K*(892)t] 
to D° — > KgT]Tr . The amplitude and phase for ao(980)Kg are fixed to 1 and 0° respectively. The 
uncertainties are statistical from the fit. "FF(Add.)" means the sum of the fit fractions for all 
modes in addition to ao(980)K^ and if* (892)77 in the model. The last row shows averaged values 
with statistical uncertainties and half the range among the four decay models. 



Additional Mode(s) 


°7Y* (892)7? 




FF(a (980)if£) 


FF(K* (892)t/) 


FF(Add.) 


NR 
^(1430)77 
#3(1430)77 + a 2 (1320)lTg 


.234±.035 
.237±.032 
.253±.031 
.269±.032 


260±10 
258±10 
251±15 
262±11 


1.350±0.097 
1.322±0.070 
1.042±0.146 
1.050±0.060 


.301±.071 
.301±.070 
.273±.050 
.310±.060 


.288±.113 
.360±.115 
.316±.097 
.186±.056 


Average and {Variation} 


.249±.032 {.018} 


259±12 {6} 


1.187±0.093 {0.154} 


.293±.062 {.019} 


.246±.092 {.087] 
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FIG. 3: The three projections of the Dalitz plot. The fit shown has contributions from ao(980)ifS, 
K*(892)r], #3(1430)77, and a 2 (1320)A"g. 

We consider possible sources of systematic uncertainties due to the background, the effi- 
ciency, the finite detector resolution, the parameterization of the matrix element amplitude, 
and choice of decay model. Central values are taken as the statistical weighted mean of the 
results summarized in Table 

For the background and the efficiency we perform the fit with the two-dimensional poly- 
nomial coefficients allowed to float constrained by their covariance matrices. We also hold 
the efficiency constant across the Dalitz plot. Deviations from the standard fit are treated 
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as systematic uncertainties. The effects are small. 

As a consistency check, we allow the parameters of one of the clearly observed resonances 
to float and extract values from the fit. For both the K*(892) and ao(980) we obtain masses, 
and width parameters consistent with previously measured values. 

Our mass resolution, small compared to the widths of the resonances we are considering, 
is a negligible effect as we observe no change when we do a fit that smears each resonance 
by a two dimensional Gaussian with widths given by propagating uncertainties on track fits 
and shower reconstructions. 

We also consider variations in the description of the decay amplitudes. We vary the 
radial parameters for the intermediate resonances between zero and twice their standard 
value of ~ 3 GeV -1 [9]. We allow the masses and widths for the intermediate resonances to 
vary within one standard deviation of their measured values |6[. The largest variation from 
the standard fit of each fit parameter is taken as an uncertainty. These uncertainties are 
combined quadratically to give a systematic uncertainty. 

The largest systematic uncertainty results from choice of decay model. Using the four 
models giving good fits we take half the range of central values, shown in Table EJ as this 
uncertainty and report it separately. 

Our analysis apparently contradicts a result done with an earlier version of our detector 
BR(D° -> K* (892)?7) = (1-8 ± 0.4)% 0, Q Q. That analysis, which focused on a search 



for this mode, made helicity angle and 77 momentum selections that are not compatible with 
a Kgr)7T° Dalitz analysis. Thus the effects of interference were not considered. Comparing 
the fit result to the Kir mass spectrum in |l4| with results obtained in this analysis we find 
that the D° — > K (892)?7 rate is larger by roughly a factor of two. 

In conclusion, we have observed for the first time the decay D° — > Kgrjir . We have 
measured the ratio of the branching fractions, 

BR(D° -> ifW) 

- g A_|_i = 0.46 ±0.07 ±0.06, (1) 

where the uncertainties are statistical and systematic respectively. Using the known D° — > 
Kgir°(K 7r ) decay rate we measure the branching fraction 

BR(D° -> K°r]TT°) = (1.05 ± 0.16 ± 0.14 ± 0.10)%, (2) 

where the final uncertainty is associated with the D° — ► K°n branching fraction. 

We have analyzed the resonant substructure of the decay D° — > Kgi]7c° using the Dalitz 
technique. We find dominant contributions from cio (980)^5 and K*(892)rj intermediate 
states. Using an isobar model including K*(892)r], do (980) lf<y, and averaging over four 
consistent models for additional components we find the amplitude, phase and fit fractions 

a K *(892) v = 0.249 ± 0.032 ± 0.013 ± 0.018, 
^(892)^ = (259 ±12±9±6)°, 
FF(K*(892)r]) = 0.293 ± 0.062 ± 0.029 ± 0.019, 

FF(a (980)K° s ) = 1.19 ± 0.09 ± 0.20 ± 0.16, (3) 

where a aoK o and f ao K° are fixed to one and zero respectively. The uncertainties are statis- 
tical, systematic, and decay model choice respectively. We also find that contributions from 
a (980)Kg and i^*(892)r7 are not sufficient to describe our data. We estimate the fit fraction 
of any additional component as 

FF(Add.) = 0.246 ± 0.092 ± 0.025 ± 0.087, (4) 
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with the uncertainties meaning as above. 
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